Biphoton generation in quadratic waveguide arrays: A classical optical simulation 
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Quantum entanglement, the non-separability of a multipartite wave function, became essential in 
understanding the non-locality of quantum mechanics. In optics, this non-locality can be demon- 
strated on impressively large length scales, as photons travel with the speed of light and interact 
only weakly with their environment. With the discovery of spontaneous parametric down- conversion 
(SPDC) in nonlinear crystals, an efficient source for entangled photon pairs, so-called biphotons, be- 
came available. It has recently been shown that SPDC can also be implemented in nonlinear arrays 
of evanescently coupled waveguides which allows the generation and the investigation of correlated 
quantum walks of such biphotons in an integrated device. Here, we analytically and experimentally 
demonstrate that the biphoton degrees of freedom are entailed in an additional spatial dimension, 
therefore the SPDC and the subsequent quantum random walk in one-dimensional (ID) arrays can 
be simulated through classical optical beam propagation in a two-dimensional (2D) photonic lattice. 
Thereby, the output intensity images directly represent the biphoton correlations and exhibit a clear 
violation of a Bell-type inequality. 

PACS numbers: 42.50.-p, 42.82.Et, 42.65.Lm 



As proven by Bell [1], nonlocal quantum correlations 
play the central role in the understanding of the famous 
Einstein- Podolsky- Rosen gedankenexperiment [2] . In the 
realm of optics, the most prominent source of such non- 
local correlations are entangled photon pairs, so-called 
biphotons [3] . They give rise to various applications such 
as quantum cryptography [4] , teleportation [5] , and quan- 
tum computation [6] . A particular robust approach to re- 
alize strong quantum correlations of path-entangled pho- 
tons in compact settings is their propagation in optical 
waveguide arrays [7, 8] that provide a unique tool for the 
experimental analysis of a spatially discrete, continuous- 
time, quantum walks with highly controllable parame- 
ters. However, in these settings the photons are typ- 
ically generated before they are launched into the ar- 
ray, usually via spontaneous parametric down-conversion 
(SPDC) in bulk optical components [3]. Recently, it was 
suggested theoretically [9] that biphotons can be gener- 
ated directly in a quadratically nonlinear waveguide ar- 
ray through SPDC, and quantum walks of the generated 
photons can give rise to non-classical correlations at the 
output of this monolithic integrated optical device. 

Remarkably, various quantum phenomena can be sim- 
ulated via purely classical light propagation in waveg- 
uide arrays, serving as an optical test-bed of fundamen- 
tal physical effects without the need of an intricate sin- 
gle photon setup [10]. In particular, the non-classical ID 
evolution of correlated quantum particles can be mapped 
onto a classical 2D evolution of a single wave packet, as 
was predicted recently by Longhi [11, 12]. By increas- 
ing the dimensionality of the structure, the same dynam- 
ics is obtained as if the number of participating photon 
wavepackets were effectively doubled. This concept was 



recently applied in the context of a discrete-time random 
walk of coherent light [13]. In this article we will show 
that linear propagation of classical light beams can be 
used to simulate the nonlinear effect of SPDC leading to 
photon pair generation and quantum walks of the gener- 
ated biphotons. Specifically, we prove theoretically and 
demonstrate experimentally that the quantum correla- 
tions of biphotons generated in a quadratically nonlinear 
ID waveguide array with pump waveguides at the edges 
can be mapped onto the linear propagation in a specially 
designed 2D system. This analogy enables us to simu- 
late, with classical light, the quantum coincidence counts 
for biphotons and the breaking of a Bell-type inequality. 



I. RESULTS 

We first extend the theory of photon-pair generation in 
quadratic waveguide arrays from Ref. [9] to finite waveg- 
uide arrays. This is especially important for an exper- 
imental realization, where the number of waveguides is 
naturally limited. We consider a type-I SPDC process 
with continuous- wave or narrow-band pump of frequency 
cop, which converts pump photons into pairs of signal and 
idler photons with frequencies uo s ^. We assume spec- 
tral filtering at half of the pump frequency, such that 
w s ,i = oo p /2. The signal and idler photons can tunnel 
between the neighboring waveguides at a rate character- 
ized by the coefficients C s ^ = C, whereas coupling of the 
pump can be neglected, due to the strong dispersion of 
the evanescent coupling [7-9, 14, 15]. Taking this into 
account and following Refs. [16-18] one can derive the 
differential equations for the biphoton wave function \I/ 
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(see Methods) 
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(1) 

whereby n s ^ label the waveguide index for signal and 
idler, respectively, S denotes the Kronecker delta func- 
tion and Af3^ is the phase mismatch in a single waveg- 
uide. The pump field is distributed over the waveguides 
according to A n whereas d e s is the effective nonlinear co- 
efficient. This equation governs the evolution of a bipho- 
ton wave function in a quadratic nonlinear waveguide 
array. The term in square brackets on the right-hand 
side describes the quantum walks, i. e. the biphoton tun- 
neling [8, 9] between the waveguides, and corresponds to 
the model derived in Ref. [12] in case of a linear waveg- 
uide array. In addition, our model also describes the 
nonlinear process of biphoton generation through SPDC, 
which is introduced by the last term. Note that both 
photons are always created at the same spatial location, 
which is expressed mathematically through (5-functions, 
but the biphotons subsequently spread out through cor- 
related quantum walks. 

After establishing the key model Eq. (1), we now 
demonstrate how the biphoton wave function can be sim- 
ulated classically. Note that Eq. (1) is formally equivalent 
to the coupled-mode equations for classical beam prop- 
agation in a square linear 2D waveguide array [12, 14] 
with additional source terms describing the pump. In 
the following, we demonstrate that for a pump exclu- 
sively coupled to the outermost waveguides of the ID 
array, the effect of these terms can be simulated clas- 
sically by introducing additional waveguides close to the 
corners of an equivalent 2D array, as shown in Fig. 1. We 
consider the case when the additional waveguide modes 
only weakly overlap with the modes at the array cor- 
ners, and the corresponding coupling coefficient C c is 
very small compared to the array, C c « C. If the light is 
launched only into the additional waveguides, then in the 
undepleted pump approximation the classical field evolu- 
tion in those waveguides is ^ieft(^) = ^ieft(O) exp(z/?i e f t z) 
and bright = bright (0)exp(z/3 righ t^), where ^i e ft(0) and 
bright (0) are the input amplitudes, f3\ e f t and /3 r i g ht are 
the mode detunings of the additional waveguides with 
respect to a single waveguide in the array. Here, we con- 
sider the case where both additional waveguides are iden- 
tical to the guides of the lattice (/3i e ft = bright = 0)- Then, 
the classical field evolution in the array follows Eq. (1) 
with the effective pump amplitudes in the edge waveg- 
uides A\ = C c ^ieft(0), An = C c bright (0) (there is zero 
effective pump, A n = 0, for 1 < n < iV), and the effective 
detuning = 0. This corresponds to perfectly phase- 

matched SPDC in the ID-array, which has been shown 
to produce the most pronounced quantum correlations 
[9]. Based on the established mathematical equivalence 
classical light propagation in the linear array structure 
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Figure 1. Settings (a) Sketch of a ID quadratic nonlin- 
ear waveguide array containing N = 7 waveguides with two 
pumps coupled to the edge waveguides leading to bipho- 
ton generation via SPDC, and the output photon detectors, 
(b) Cross-section of a 2D waveguide array for the classical 
optical simulation of biphoton generation via SPDC consist- 
ing of 3 x 3 waveguides. The additional waveguides represent 
the pump beam coupled to the edge waveguides. 



shown in Fig. 1(b) can be employed to simulate the quan- 
tum properties of biphotons generated in a quadratic 
nonlinear array shown in Fig. 1(a). In particular the 
photon number correlation T n ^ ni = («J ls aJ l .a n .a ns )^, 
which describes the probability of simultaneous detec- 
tion of photons at the output of the waveguides with 
number n s and is of interest. Thereby, a ns(i) denotes 
the photon annihilation operator in these waveguides. In 
the classical simulation the output intensity distribution 



corresponds directly to T: 



F( cl ) - T 

L n s ,m ~ 1 n s ,n i 



\n s ,n ? 



In s ,rii 



^n sy rii • (2) 



To verify the correct operation of our simulator, we 
first implemented a system with a single outer pump 
waveguide coupled to the corner of 3 x 3 waveguide array 
[Fig. 2(a)]. For the investigation light is launched into 
the device and imaged onto a CCD at the output (see 
Methods). The experimental output intensity distribu- 
tion is shown in Fig. 2(b). We perform a simple scaling 
of the measured intensity according to Eq. (2) and de- 
termine the two-photon correlation function as shown in 
Fig. 2(c). These results closely match the theoretical pre- 
dictions presented in Fig. 2(d). 

After verifying the correct operation in the single- 
pump regime, we return to the initially considered case 
of two pump waveguides at the corners of the 2D lattice 
as shown in Fig. 3(a). The challenge here was to have 
light of the same amplitude and phase in both pump 
waveguides. The implementation is discussed in Meth- 
ods. We likewise use the output intensity distribution 
[Fig. 3(b)] to calculate the photon number correlation 
[Fig. 3(c)] which agrees very well with the simulation 
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Figure 2. Results for one pump waveguide (a) Sketch 
of the linear 2D waveguide array simulating quantum corre- 
lations of biphotons generated by SPDC with one pump at 
the edge waveguide, (b) CCD camera image of the output 
light distribution when only the pump waveguide is excited. 
(c,d) Correlation map of the simulated ID quantum system: 

(c) Extracted from intensity measurement shown in (b) and 

(d) calculated numerically based on Eq. (1). (e,h) Simulated 
non-classicality function determined with Eq. (3) for (e) ex- 
perimental and (f) numerical correlations from (c) and (d), 
respectively. 



shown in Fig. 3(d). 



II. DISCUSSION 

After obtaining the correlation functions from 
the output intensities, we determine in ac- 
cordance with Ref. [8], the similarity S = 



; V/( E 

the numerically calculated (T, 
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obtained (r^ p 



num\ 

distributions as S 



3 i,3 

and experimentally 
= 0.954 for the case 
of one pump waveguide and S = 0.972 for two pump 
waveguides. 

In a further step we demonstrate that our classical op- 
tical simulator can model non-classical biphoton statis- 
tics. Following Ref. [8], we calculate the non-classicality 



\/"^ n s ^Tii ,rii ^n s 



(3) 



Values of V n ^ ni > 0, indicate true quantum behavior cor- 
responding to a violation of Bell's inequality, which can- 
not occur in a purely classical setting. In the single pump 



Figure 3. Results for two pump waveguides (a) Sketch 
of the linear 2D waveguide array simulating quantum corre- 
lations of biphotons generated by SPDC with equal pumps at 
the two edge waveguides, (b-f) Experimental and theoretical 
results, notations are the same as in Fig. 2. 



case we plot the function V nsini for our experimental data 
in Fig. 2(e), which shows the presence of positive values 
in agreement with numerical simulations [Fig. 2(f)]. This 
proves that our 2D classical setting successfully simulates 
actual non-classical ID statistics. Fig. 3(e) and Fig. 3(f) 
illustrate the experimentally and numerically calculated 
non-classicality. Also for this setting it is evident, that 
the nonclassical features can be simulated with high qual- 
ity. 

As evident from the strong overexposure of the CCD 
at the pump waveguide in Fig. 2(b) and fact that the 
majority of the light remains in the pump waveguides in 
Fig. 3(b), the undepleted pump approximation holds. 

For both a single and double pump waveguides, the 
output correlation should be symmetric, T nsjfl . = r n . >ns , 
due to the indistinguishability of two photons, and we 
see that this feature is properly reproduced through the 
2D experimental intensity distributions where I na ,m — 
Ini,n s i see Figs. 2(b) and 3(b). For a double pump config- 
uration, there appears additional symmetry due to identi- 
cal pump amplitudes at the two boundaries, and we have 
in theory I nsjfli = lN+i- nij N+i-n s [Figs. 3(d)], whereas 
such symmetry is clearly broken for a single pump case 
[Figs. 3(d)]. In experiment we have only a slight asym- 
metry of the output intensity distributions due to imper- 
fections of the three-waveguide input-coupler (see Meth- 
ods). Nevertheless, for all cases we have similarity S close 
to unity, as well as a good agreement of the experimen- 
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tal and simulated non-classicality which proves excellent 
accuracy of the classical optical simulator. 

In conclusion, we have derived a model describing the 
evolution of a biphoton wave function in ID quadratic 
nonlinear waveguide arrays, where the photons are gen- 
erated through SPDC and undergo correlated quantum 
walks. We further demonstrated analytically as well as 
experimentally that the quantum biphoton dynamics can 
be simulated by a classical wave evolution in a linear 2D 
waveguide array with additional waveguides represent- 
ing the pump. Hereby, the additional spatial dimension 
provides the degrees of freedom which are otherwise en- 
coded in the two-particle dynamics. The classical mea- 
surements of the output light intensity directly simulate 
the quantum biphoton correlation function, in particu- 
lar including the regime of the breaking of a ID Bell-like 
inequality. 



III. METHODS 

A. Derivation of the biphoton wave function 

In the following a detailed derivation of the differential 
equation for the biphoton wave function is presented. We 
start from the coupled mode equations for the classical 
light amplitudes in waveguide n (1 < n < N) 



i^ + C(E n _ 1+ E n+1 ) 



0, 



(4) 



with the boundary conditions E$ = E^+i = 0. The 
eigenstates of the system are the supermodes En(z) = 
exp(i/3 m z), where m is the mode number (1 < 
m < N), j3 m = 2Ccos(/c m ) is the propagation constant, 
= y/2/(N + 1) sm(nk m ) is the amplitude distribu- 
tion, and km — irm/(N + 1). For the pump, the cor- 
responding coupling coefficient C p would generally have 
a much smaller value compared to the signal and idler 
waves, C p <C C, due to the weaker mode overlap be- 
tween neighboring waveguides at higher frequencies [14]. 
Therefore, we neglect coupling effects for the pump beam 
(C p « 0), and assume its amplitude A n to be constant 
along the propagation. 

Subsequently, we employ the mathematical approach 
of Refs. [16-18] developed for multi-mode waveguides, 
and obtain the following expression for the two photon 
state |S): 

N N 
rrii — l m s = l 



where 



^m s ,mi(z) =ex.p(if3 ms z + if3 mi z) 



[ d^7 ms , m . exp(iA/3 ms , m V). 
Jo 



(5) 



n ( 6 ) 



Here B is a constant, a) are photon creation opera- 
tors at the numbered supermode states, |0,0) is the 
vacuum state, \I/ ms , mi is a two-photon wave function, 
A/3 ms?m . = A/^K ) — j3 rris — f3 mi is a phase mismatch be- 
tween the modes, and is the mismatch in a single 

waveguide. The value 7 ms , m . = Eli 4ff A^i™ 5 ^™^ is 
the normalized spatial overlap of the signal and idler dis- 
tributions with the pump distribution. We note that an 
integration over frequency is omitted, since we consider 
a frequency window narrower then the phase-matching 
bandwidth of the SPDC. 

For our analysis, it is convenient to rewrite Eq. (6) in 
the form of differential equations, 



1 dz 



■-(Pm a + Pm i )^m s ,m i 
+ 11m s ,m i exp(z 



(7) 



This allows us to change our representation to the real- 
space wave function, which is related to the modal formu- 
lation as ^n s ,n, = Em^m, ^m^m^T^!?^ • D ^e to the 

orthogonality of the supermodes, we can invert these re- 
lations as ^m s ,mi = J2n s ,m ^ n s .n^^ . Substitut- 
ing the latter expression in Eq. (7) and again employing 
the orthogonality property together with the fact that 
by definition the supermodes satisfy Eq. (4), yields the 
differential equations (1) for the biphoton wave function 
in real-space representation. 



B. Experimental realization 

To realize the classical light simulator of quantum 
biphoton statistics, we fabricate waveguide arrays in 
fused silica using the femtosecond direct-write technique 
[19, 20] where we used 160 fs pulses with an average power 
of 20 mW and a writing velocity of 60 mm/min. A sketch 
of the writing setup is shown in Fig. 4(a). Our optical 
simulator consists of a 3 x 3 waveguide array with outer 
pump waveguides coupled to the corners of the array as 
shown in Fig. 4(b). To ensure homogenous horizontal 
and vertical coupling we rotated the structure by 45° 
and fabricated the waveguides in a rhombic geometry. 
To meet the undepleted pump approximation of SPDC 
we made sure that the intensities in the pump waveguide 
do not decrease significantly by choosing a sufficiently 
weak coupling of C p = 0.125 cm -1 , corresponding to a 
distance between the pump waveguide and the array of 
25 fim [15]. Laser light of 633 nm wavelength was in- 
jected into the pump waveguides and the output inten- 
sity (In s ,ni) was observed with a CCD camera. We then 
use the measured intensity to determine the simulated 
biphoton correlation. 

In a first step we investigated a simulator having a sin- 
gle pump waveguide [Fig. 2(a)]. The distance between 
the array waveguides is 17 /im corresponding to a cou- 
pling constant of C = 0.78 cm -1 . 
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Figure 4. Fabrication setting (a) Sketch of the writing 
setup, where femtosecond laser pulses are focused into a trans- 
parent bulk material. In the focal region, the refractive index 
of the material is permanently increased, (b) Sketch of the 
setting with two edge pump waveguides, where the phase- 
matched symmetric amplitude distribution of the pump is 
achieved by a specially designed three-waveguide directional 
coupler. The light is launched into the central waveguide 
marked with an arrow. Additionally, a micrograph of the 
output facet is shown as inset (scale bar equals 40 /im). 
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